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1^ • The masses of 1 charmonium and bottomonium hybrids 

are evaluated in terms of QCD sum rules. We find that the 
ground state hybrid in charm sector lies in mn^ = 4.12 ~ 4.79 
^ I GeV, while in bottom sector the hybrid may situated in rriH^ = 

10.24 ~ 11.15 GeV. Since the numerical result on charmonium 
'nI" ' hybrid mass is not compatible with the charmonium spectra, 

■ including structures newly observed in experiment, we tempt 

to conclude that such a hybrid does not purely exist, but 
^ ! rather as an admixture with other states, like glueball and 



regular quarkonium, in experimental observation. However, 
our result on bottomonium hybrid coincide with the "exotic 
^ ■ structure" recently observed at BELLE. 



ix^ ■ 1 Introduction 



Quantum chromodynamics (QCD) is believed to be the underlying theory for strong 
interaction. It is commonly believed that the mechanism responsible for the hadronic 
properties is subject to the non-perturbative aspect of QCD. Unlike the perturbative 
QCD which is well understood in some sense, we do not have a reliable and effective way 
to tackle with the non-perturbative QCD effect yet. In this respect, to get a deep insight 
in the physics associated with the non-perturbative QCD is one of the most important 
tasks for the society of high energy physics. So far and in this aim, to evaluate the physical 
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quantities of hadrons, such as hadron spectra, hadronic transition matrix elements, parton 
distributions and fragmentation functions, in most cases people have to invoke to typical 
phenomenology models. 

Apart from the so-called regular hadrons, the meson and baryon, QCD also does 
not exclude exotic hadronic structures, like hybrid, glueball and multiquark structures. 
Normally, the hybrid meson refers to a state which contains a pair of constituent quarks 
and a dynamic gluon. In the color-fiux tube model, instead, the hybrid corresponds to the 
structure where the gluonic degree of freedom is excited. Even though these two pictures 
look different, they may be just two sides of the same object. Thus, whether they are 
reconcilable with each other or not is an interesting question. Indeed, a careful study 
may shed light on the hybrid structure and help us to get a better understanding of the 
non-perturbative QCD effects. 

Up to now, a many of effective methods was proposed [11 [2], [Sj H] in evaluating the 
hybrid mass spectrum, for instance, by studying the quarkonium hadronic transition via 
multipole expansion [5l[6], the bag model [7], the flux-tube model [8j, lattice QCD |9] and 
QCD Sum Rules [IOl[IIl[l2l[l3l|l2[T5l[T6l[I7l|18]. Whereas, those theoretical evaluation 
results diverse greatly with each other, and hence it is hard to pin down any exotic 
structures as hybrids in experiment. Therefore, further theoretical investigations are 
necessary, and they were partially done. For instance, the light hybrid masses evaluated 
in [19] were updated by Narison and his collaborators [201 l2T| [22] in the framework of 
QCD Sum Rules. Moreover, as indicated in Refs. [23l [2l] . the hybrid states may not 
exist independently, but may rather admixtures of hybrids with regular quarkonia or 
even glueballs with the same quantum numbers. To clarify the messy picture, a more 
accurate evaluation of pure hybrid spectrum is still vital, namely one can then confront 
the theoretical predictions with the experimental data to determine the hybrid component 
of a hadronic state. 

Among those theoretical methods in dealing with the non-perturbative effects, QCD 
Sum Rules innovated by Shifman et a/.[2S] turns out to be a remarkably successful and 
powerful technique for the computation of hadronic properties. By virtue of QCD Sum 
Rules, hybrids with various quantum numbers and the flavors have been investigated. 
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For light-quark hybrids, in order to avoid the mixing between hybrids and ordinary 
mesons, Ref.[TO] considered specifically the hybrids possessing exotic quantum numbers 
1 '", and Ref.[TT] took another exotic quantum numbers into consideration, and ob- 
tained the relative masses and decay amplitudes. Employing the heavy quark effective 
theory (HQET), Huang, Jin and Zhang evaluated the masses and decay widths of several 
typical heavy- light hybrids [15]. The heavy-quark hybrids masses were also evaluated in 
Refs. [161 HH [18] through QCD Sum Rules. 

Recently, a hadronic structure with mass 4324 ±24 MeV and width 172 ±33 MeV, has 
been observed by the BABAR Collaboration in the iIj{2S)7t^'k^ invariant-mass spectrum 
[26] . This structure is obviously different from y(4260) reported in Ref.[27]. Later on, two 
resonant-like structures are observed in also the iP{2S)tt^tt~ mode by the Belle Collabo- 
ration, one resides at 4361 ±9 ±9 MeV with a width of 74± 15 ± 10 MeV, which coincides 
with what reported by the BABAR collaboration, and another is at 4664 ± 11 ± 5 MeV 
with a width of 48 ± 15 ± 3 MeV [28J. A variety of theoretical speculations has been put 
forward for these hadronic structures. For instance, supposing Y(4664) still be a normal 
member of the charmonium family, it was interpreted as different charmonium states, like 
a 5^ Si state, a 6^ Si state, or even a 5^5"! — i^Di mixed state [29]. The F(4664) was also 
interpreted as a baryonium state [30], the radial excited state of -^(lAcAc) ± 
and a /o(980)^/'' molecule [31j. Starting from the QCD Sum Rules, Albuquerque et al. 
[32] computed the mass of F(4664) based upon the assumption that it is a vector cscs 
tetraquark state. In the literature, for convenience, y(4361) and y(4664) are usually 
named as y(4360) and y(4660), respectively. 

Since a series of newly observed "exotic" states in charmonium energy region is 1 
hadrons, in this paper we reinvestigate the 1 charmonium hybrid, which is constructed 
by a pair of charm-anticharm quarks and a gluon, by means of the QCD Sum Rules. 
In our calculation, the interpolating current of the charmonium hybrid is chosen to be 
Qsip'j'^'j^T^Gtplx), which can be easily found having the correct quantum numbers of the 
hybrid and was also used in Ref . [TB] , where, however, the tri-gluon condensate contribution 
was not taken into account and the plane wave method was used. In our work, we keep 
the operator product expansion (OPE) to dimension six, the dimension of the tri-gluon 
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condensate, and take the widespreadly used Fixed-Point gauge technique. Our numerical 
result indicates that the contribution of the tri-gluon condensate is not negligible and 
even somehow important for the estimation of the charmonium hybrid mass. In a similar 
work done by Kisslinger, Parno, and Riordan [18], the current J^i^) = ^C'~fi,G'^'^'^{x) 
was employed to calculate the 1 charmonium hybrid, which we think is improper due 
to the incompatible quantum number with the concerned hybrid, and they obtained a 
quite low mass of 3.66 GeV. 

The rest of this paper is organized as follows. In Sec. II we derive the formulas of the 
correlation function U^^, in terms of the QCD Sum Rules with the interpolating current 
for J^^ = 1 . In Sec III, our numerical evaluations and relevant figures are presented. 
Section V is remained to summary and concluding remarks. 



2 Formalism 



In the framework of QCD Sum Rules, the starting point is properly constructing the 
correlation function, i.e., 

= ^ J d'xe''^-^{0\T{J^ix)Jlm\Q) . (1) 

Here, the interpolating current for the heavy hybrid with quantum number J^"-^ = 1 
is chosen to be 

Ux) = g,rix)Yi,^G;,{x)^\x) , (2) 

where Qs is the strong coupling constant, a and b are color indices. A" is the color matrices, 
and G'^^{x) = e^,yaisG^'"'^{x)/2 is the dual field strength of G'^j^{x). Generally, the two- 
point function 11^;^ may contain two distinct parts, the vector part ny(g^) and the scalar 
part Usiq"^) which represent the contributions of the correlation function to the vector 
channel J^'" = 1 and scalar channel J^*" = 0+^, respectively. They can be explicitly 
expressed as 

n,.(g) = {^-^-9,.)Ily{q') + ^-^Us{q'). (3) 
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Since the main task of our work is to study the mass of the vector heavy hybrid with the 
quantum number 1 , following we only analyze the vector part Ilv{q^)- 




c D 



Figure 1: The typical Feynman diagrams for the calculation of heavy vector hybrid masss. 
Here, the permutation diagrams are implied. Diagram (A) represents for the contribution 
of unit operator; (B) for the contribution of two-gluon condensate; (C) and (D) for the 
contribution of three-gluon condensate. 



By the operator product expansion (OPE), the correlation function ny(g^) can be 
written as 

Uv{q')^U^''\q') + Ur''{q'), (4) 

where to satisfy the necessity of our calculation, W^^{q'^) is obtained by taking the imag- 
inary part of the Feynman diagram A, and 11^°^'^ (q^) represents the contributions coming 
from all possible condensates. In this work, we consider the condensates up to dimension 
six. 

First, we calculate the imaginary part, the absorptive part, of the Feynman diagrams 
which represents the perturbative contribution to the correlator, and the result reads 

,6 



asm} 



15r + 185r - 778^'^ - 496r + 12m - 192 



+15t^Vl^t{t^ - I2f + m - 128) log ^ (5^) 

y/t -I 

where, t = imq/s, and mg is the mass of the heavy quark, and p^^^^{t) = Imn(t). 
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The contributions of non-perturbative condensates from diagrams in Figure 1 are 
nr''^(g') = ^'dx^^{[8(l-x)V-llm|]+ln(A)[2(l-x)V-3m|]}, (6a) 
nr^-(0 = iyj-^[3xHA) + '-^ + 17x], (6b) 



(^?'^')ro o M [2(3- l,r)///5 + ,r(ll,r2-27,r + 13)^/2]x 



A 



ncond,D(^2) ^ 2 da;^gl{2a;(2-3a;)ln(A)- 
(x - I)g2[3xg2(a; - 1)2 + (2 - 3x)m|]x2 2(5-24x)x, , , 
+ + 3 }■ (6c) 

Here, A = —{l — x)xq'^ + mQ, and symbols B, C, D represent the corresponding diagrams, 
respectively. 

In order to eliminate contributions from higher excited and continuum states, a well- 
known procedure, the Borel transformation, is performed to above obtained results, and 
then we get 

B[nr''^(g2)]) ^ du dx j' ' ^ . 2e"^^^^ - S-^e , (7b) 

Jo Jo 1927r^(l -x)\ lo^ ) 

+e ^b(i--)-(13x- ll)a;^] +e ^b(i-)- (6x - 5)a;^m| ^ , (7c) 



where, Mb is the Borel parameter. 

Suppose the existence of the quark-hadron duality, the resultant sum rule for the mass 
of the vector heavy hybrid reads 



with 



^H^\l-^ (8) 



Ro^- r pP^'\s)e-'/^B + B(nr'^'^)+B(n^6°"^'^) + B(n^°"^'°) , (9) 

d 

Ri = — ^Ro- (10) 



Here, sq is the threshold cutoff introduced to remove the contribution of the higher excited 
and continuum states [331. 



3 Numerical Analysis 



For numerical calculation, the leading order strong coupling constant 

Alt 

c.(Ml) = — (11) 

is adopted with Aqcd = 0.220 GeV and nj being the number of active quarks. Of the 
two- and three-gluon condensates we take both the prevailing values 



(a^G^) = 0.038 ± 0.011 GeV^ , {g^,G^) = 0.045 GeV^ (12) 

and also recently obtained ones [M], i.e., 

{asG^) = 0.075 ± 0.020 GeV^ , {g^G^) = (8.3 ± 1.0)GeV2 (a.G^) (13) 

into account. The heavy quark masses are taken to be: 

TJic = 1.26 ~ 1.47 GeV, nib = 4.22 ~ 4.72 GeV . (14) 

Here, the masses span from the running masses in MS scheme to the on-shell masses of 
QCD Sum Rules 



For the selection of an appropriate Borel parameter M^, we adopt the criteria proposed 
in Refs.[2Sl ESI EZ]- Defining mjy(M|) in Eq.§i) to be fthcovv{M^) while the continuum 
contribution is absent, i.e., (sq = oo), and mniM'^) to be mH,nopower(^B) in case of no 
power corrections. Then we get two discrimination functions /cont(^|() and /nopower(^B)5 
satisfying 

/c„„.(M|) = P^,, (15) 

'"'H,nopowcrl,-'^^ B) 
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According to the sum rule criteria, the mass function obtained in Eq.([8]) is vahd only in 
the situation of being neither too small nor too large. In case is very small, the 
omitted terms of high- dimensional condensates, which are proportional to high powers of 
may become too important to be neglected. To get a reliable prediction of the 
hybrid mass in QCD Sum Rules, /nopowerl^l) is required to be less than 10% from unit, 
which ensures the contributions from the non-pertubative condensates to be much less 
than what from the perturbative term. On the other hand, a very large implies the 
invalidation of quark-hadron duality approximation. Normally, the /cont(^B) is required 
to be 70% more to suppress the contributions of higher resonances and continuum. 



A B 




2468 10 2468 10 



Figure 2: The curves of /nopower and /cont for charmonium hybrid versus Borel parameter 
while the continuum threshold cutoff sq = 26 GeV^. Figures A and B correspond to 
the gluon condensates from [33] and [3l], i.e. Eqs.(|T2l) and ( TT3l) . respectively. 

For charmonium hybrid, to find the reliable sum rule for hybrid mass according to the 
aforementioned criteria, in Fig. [21 we draw the curves of /nopower and /cent as functions 
of M^. The figure indicates that /nopower and /cont niay satisfy the above mentioned 
requirements, i.e., the contribution from pole terms is more than 70% and the contribution 
from condensates is less than 10%, while lies in between 1.80 to 7.80 GeV^ for Eq. f[T2]) 
and between 2.10 to 7.20 GeV^ for Eq.( TT^ . respectively. 

For bottomnium hybrid, as shown in figure 3 the reliable range for lies in between 
3.00 to 34.00 GeV^ for Eq.dH]) and between 5.00 to 33.00 GeV^ for Eq.(II3]), respectively. 
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A 



B 




Figure 3: The curves of /nopower and /cont for bottomnium hybrid versus Borel parameter 
while the continuum threshold cutoff Sq = 130 GeV^. Figures A and B correspond to 
the gluon condensates from [33| and p4j|, i.e. Eqs. (fT2|) and ( IT3|) . respectively. 



With the above requirements in establishing the sum rule, to calculate the physical 
quantities in terms of QCD Sum Rules, one needs to find an optimal window for Borel 
parameter and threshold parameter sq. Within this window, the physical quantities, 
here the hybrid mass, are maximally independent of the Borel parameter. 

In figure HJ- A, we draw lines for different threshold parameter sq, that is 29, 26, and 
23 GeV^ while Borel parameter varying from 1.80 to 7.80 GeV^ The fi gure evidently 
indicates that there exists a window for the Borel parameter between 4.20 GeV^ and 
7.80 GeV^, in which the evaluated hybrid mass is mostly independent of M^, especially 
m case of So = 26 GeV^ Namely, it is proper to select the threshold parameter to be 
26 GeV^, which hints that the mass of the first excited state is above 5.10 GeV. In figureH]- 
B, the corresponding parameters are sq = 29,26,23 GeV^ 2.50GeV2 < M| < 7.20GeV^ 
and the proper threshold parameter is also 26 GeV^. 

Considering the uncertainties remain in the input parameters, the quark mass, conden- 
sates, the Borel parameter and continuum threshold sq, we obtain the charmonium 
hybrid mass to be 

TUH^ = A.52t'ofs GeV . (17) 
Here, the charm quark mass rric goes from 1.26 ~ 1.47 GeV; the condensates take the 
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Figure 4: The dependence of 1 charmonium hybrid sum-rule mass on the Borel pa- 
rameter Mg with different continuum threshold cutoffs sq, that is 29 GeV^, 26 GeV^ 
and 23 GeV^ from up to down in the figure. Figures A and B correspond to the gluon 
condensates from [33j and i.e. Eqs. f|T2|) and f lT3|) . respectively. 



magnitudes of Eq. fll2p ; the Borel parameter varies from 4.20 GeV^ to 7.80 GeV^; and 
the continuum threshold Sq changes from 23 GeV^ to 29 GeV^. The central value of mn^ 
in Eq. (IT7|) is reached by taking the central values of the quark mass and condensates, 
while setting sq = 26 GeV^ and M| = 6.00 GeV^ 

For the gluon condensates taken as Eq. fll3l) . we obtain the charmonium hybrid mass 
to be niH, = 4.45+°:| GeV with sq = 26GeV2 and = 6.00GeV^ which is slightly 
higher than fll7p . 

In bottomnium sector, by the same procedure, but with inputs Sq = 130GeV^ and 
M| = 20GeV^ as shown in fi gure 5, we readily obtain the 1 bottomonium hybrid 
mass, that is rjiH^ = ^^-^^-oll GeV and 10.70;l^Q;4gGeV for condensates from [53] and 
[Sl|(Eqs.([I2) and (dnD), respectively. 

4 Summary and Conclusions 

In this work we recalculate the 1 heavy quarkonium masses in the framework of 
QCD Sum Rules. The central part of this calculation relys on the evaluation of the 
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Figure 5: The dependence of 1 bottomnium hybrid sum-rule mass on the Borel pa- 
rameter with different continuum threshold cutoffs Sq, that is 140 GeV^, 130 GeV^ 
and 120 GeV^ from up to down in the figure. Figures A and B correspond to the gluon 
condensates from j33] and [M], i.e. Eqs. (fT^ and ( 1T5]1 . respectively. 



Wilson coefficients of the operators for the two point correlation function constructed with 
a suitable hybrid current. In former studies [HI [18], people also fought the same target 
but with some differences from this work. In previous works, operator expansion only up 
to dimension four, but in our study the dimension six operators are taken into account. 
And our results indicate that the dimension six operator, the three-gluon condensate, is 
important in attaining a wide stable plateau and hence a stable sum rule, which makes 
the predictions for the masses of ccG and bbG hybrid states more reliable. We find 
the interpolating current used in Ref.[18j is improper, and our procedure in establishing 
the sum rule differs from what performed in Ref.[T6]. In our calculation, the central 
values of charmonium- and bottomonium-hybrid masses are 4.52 GeV and 10.81 GeV for 
condensates in (fT2|) and 4.45 GeV and 10.70 GeV for condensates in (fT3|) . respectively. 

Considering of the Y states in charmonium region which are observed recently at the 
B-factories, since our predicted mass of the charmonium hybrid resides between y(4360) 
and y(4660), we are tempted to conclude that neither of the y(4260), y(4360) and 
y(4660) states can attribute to pure charmonium hybrid state. If we take the errors into 
consideration seriously, the masses of y(4360) and y(4660) are closer to our estimation, 
thus might be candidates. However, at present stage it is still hard to draw a definite 
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conclusion. As a matter of fact, by the discussion given in the introduction, a pure hybrid 
or glueball might not independently exist. The observed structures could be admixtures 
of relevant hybrid with glueballs and regular quarkonia with the same quantum numbers 
[23l 121]. Therefore, in this sense it is understandable why the theoretical hybrid mass 
does not coincide with any of the resonances observed in experiments. To clarify this 
issue, i.e., to calculate the mixing angles among states with the same quantum numbers, 
we need a bigger database on the exotic states, which might be available at the LHCb 
and planned Super-Flavor factory. 

Finally, we have also estimated the mass of 1 bottomonium hybrid state with quite 
small uncertainties, the bbG, which is compatible with the recent BELLE observation of 
the exotic y;,(10890) [38], and may confront to the LHCb data in near future. 
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